The eccentric universe 



Arjun Berera* , Roman V. Buniyf and 
Thomas W. Kephart| 

* School of Physics, University of Edinburgh, Edinburgh, EH9 3JZ, United 
Kingdom 

t Institute of Theoretical Science, University of Oregon, Eugene, OR 97403, USA 
X Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 
37235, USA 

E-mail: ab@ph.ed.ac.uk, romanOuoregon.edu and 
thomas.w.kephart@vanderbilt.edu 

Abstract. For a universe containing cosmological constant together with 
uniform arrangements of magnetic fields, strings, or domain walls, exact solutions 
to the Einstein equations are shown to lead to a universe with ellipsoidal 
expansion. The magnetic field case is the easiest to motivate and has the highest 
possibility of finding application in observational cosmology. 
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1. Introduction 

The first year WMAP results [1, 2, 3] contain interesting large-scale features which 
warrant further attention [4, 5]. First, there is the suppression of power at large angular 
scales {6 > 60°), which is reflected most distinctly in the reduction of the quadrupole 
C2. This effect was also seen in the COBE results [6, 7]. However, WMAP was able 
to make more precise measurements of C2. Next, the quadrupole C2 and octupole 
C3 are found to be aligned. In particular, the £ = 2 and 3 powers are found to be 
concentrated in a plane inclined about 30° to the Galactic plane. In a coordinate 
system in which the equator is in this plane, the £ = 2 and 3 powers are primarily in 
the m — ±£ modes. The axis of this system defines a unique direction and supports 
an idea of power in one direction being suppressed compared to the power in the 
orthogonal plane. These effects seem to suggest one (longitudinal) direction may have 
expanded differently from the other two (transverse) directions, where the transverse 
directions describe the equatorial plane mentioned above. 

After the COBE experiment, Monte Carlo studies were used [8] to cast doubt on 
quadrupole suppression and analyses of the WMAP data [1, 2, 3, 4, 5] have arrived 
at similar conclusions. Nevertheless, interesting physical effects are not ruled out and 
a realistic physical model that could explain some or all effects in the WMAP data 
would be of interest. 

It would be most satisfying to explain global anisotropy in the universe by a simple 
modification of the conventional Friedman-Robertson- Walker (FRW) model. As a first 
step in this direction, we will consider an energy-momentum tensor which is spatially 
non-spherical or spontaneously becomes non-spherical. Such situations could occur 
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when defects or magnetic fields are present. Magnetic fields [9] and cosmic defects [10] 
can arise in various ways. Moreover, large scale magnetic fields exist in the universe, 
perhaps up to cosmological scales [9, 11]. These considerations motivate us to examine 
the effect magnetic fields and defects have on the expansion of the universe [12]. 

To begin to understand the form, significance and implications of an asymmetric 
universe, we here modify the standard spherically symmetric FRW cosmology to a 
form with only planar symmetry [13]. Our choice of energy- momentum will result in 
non-spherical expansion from a spherical symmetric initial configuration: an initial 
co-moving sphere will evolve into an ellipsoid that can be either prolate or oblate 
depending on the choice of matter content. We then explore the effects of this type of 
asymmetric expansion on modes of a multipole expansion. For the sake of clarity, we 
first give some general properties of cosmologies with planar symmetry (The universe 
looks the same from all points but they all have a preferred direction.) and then 
investigate one case in detail. Our exemplar will be a universe filled with uniform 
magnetic fields and cosmological constant. This is perhaps the most easily motivated, 
analytically solvable case to consider and it will give us a context in which to couch the 
discussion of other examples with planar symmetry and cases where planar symmetry 
is broken. 

To set the stage, consider an early epoch in the universe at the onset of cosmic 
inflation, where strong magnetic flelds have been produced in a phase transition 
[14, 15, 16, 17, 18]. The magnitude of the magnetic field and vacuum energy 
(A) densities initially are about the same, but eventually A dominates. It was 
estimated [16] that the initial magnetic field energy produced in the electroweak 
phase transition was within an order of magnitude of the critical density. Other 
phase transitions may have even higher initial field values [17, 18]. Hence, it is not 
unphysical to consider a universe with B-fields and A of comparable magnitudes. If 
the B-fields are aligned in domains, then some degree of infiation is sufficient to push 
all but one domain outside the horizon. (Below we also discuss cases where there are 
still a few domains within the horizon.) 

2. Planeir symmetric cosmology 

To make the simplest directionally anisotropic universe, we modify spherical symmetry 
of FRW space-time into planar symmetry. (Cylindrical symmetry is, of course, not 
appropriate since it introduces a preferred location of the axis of symmetry.) The 
most general form of planar-symmetric metric (up to a conformal transformation) 
is [13] {gfi„) = diag(l, — e^", — e^", — e^^), where a and b are functions of t and z; and 
the xy-plane is the plane of symmetry. We also impose translational symmetry along 
the z-axis; the functions a and b now depend only on t. For the above metric, the 
nonzero-components of the Ricci tensor then are 

R°o = -{2a + b + 2a^ + b^), 
iJ^i = = -(a + + ab), 
= -(b + b'^ + 2ab). 

To support the symmetry of space-time, the energy momentum tensor for the 
matter has to have the same symmetry: 

(r%) = (87rG)-Miag(^,r7,r;,C). (1) 
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Here the energy density ^, transverse r] and longitudinal C tension densities are 
functions of time only. Examples of plane-symmetric spaces include space uniformly 
filled with either magnetic fields, static aligned strings, or static stacked walls; i.e., 
the defects are at rest with respect to the cosmic background frame. Of course, we 
can add to these any spherically-symmetric contributions: vacuum energy, matter, 
radiation (see Table 1). The entries in Table 1 are scaled according to Eq. (1); for 
instance, A = SttGA. 



Table 1. Various contributions to the components of the energy-momentum 
tensor in Eq. (1). 
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The Einstein equations corresponding to Eq. (1) are 

+ 2ab = ^, (2) 

a + b + a'^ + ab + b'^ = r], (3) 

2a + = C- (4) 

In addition, we need the equation expressing covariant conservation of the energy- 
momentum [a direct consequence of Eqs. (2)-(4)]: 

^ + 2d(e-r?) + 6(^-C) = 0. (5) 

Before considering specific models for the energy-momentum, we first establish 

several general features of an anisotropic universe described by Eqs. (2)-(5). 
We assume that before anisotropic effects became important, the universe was 
isotropic [19], aj = 6j, and was expanding isotropically, dj = 6j > 0. Without loss of 
generality we set aj = 0, which is equivalent to setting the initial scale factor to unity. 
During the initial phase of anisotropic expansion when non-uniform contributions 
to the energy- momentum are more significant, different tension densities in the 
longitudinal and transverse directions cause the universe to expand non-spherically. 
At a later time, when all contributions except for the vacuum energy fade away, 
longitudinal and transverse expansion rates become equal and the expansion again 
proceeds isotropically. Thus, an initial co-moving sphere evolves into an ellipsoid; 
whether it is oblate (a > b) or prolate (a < b) depends on which tension dominated 
during the initial phase of deformation. 

For all known forms of classical matter, the components of the energy-momentum 
tensor satisfy the dominant energy condition [20]; in our case it says ^ > 0, ^ > ry, 
and ^ > C (see Table 1 for examples). These conditions lead to various bounds and 
asymptotics of which we list here a few important examples. 

First, one can show that for the initial conditions given above, the energy density 
never increases. Also, the space always expands transversally, and unless the magnetic 
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field contribution dominates, the space expands longitudinally. Next, the transverse 
expansion rate has its maximum a,t t = and the same is true for the longitudinal 
expansion rate unless the wall contribution dominates. 

Second, the eccentricity [21] is bounded. Indeed, Eqs. (2) and (4) lead to 
a — b = (C~C~ 2d) /2a, which upon using the energy condition ^ > C gives the 
bound e""'' < dj/a. Solving Eq. (4) asymptotically for large t, we find a ~ (C/3)2, 
which leads to the asymptotic bound 

e"-'<(ei/C)^. (6) 
Note that in an expanding universe > ^ > C and so (6) does not forbid a > b. 
A sub-case with vacuum energy and any combination of magnetic fields and strings 
aligned in the same direction has ^ = C) and similarly we find 

e"-'~(ei/C)^. (7) 
Third, let us expand Eqs. (5) for large t to find 

(5C)-+n(iA)^5e = 0; (8) 

here n = 3 — 2Sr]/6^ — S(/6^. For cosmological constant plus magnetic field, strings, or 
walls, n = 4, 2 or 1, respectively. Thus the energy density approaches its asymptotic 
value A according to 6^ cx e~*/*« with characteristic time = n^^{3/X)^ . 

We consider here three exactly solvable cases: cosmological constant plus either 
uniform magnetic field, homogeneously aligned strings, or homogeneously stacked 
walls. We concentrate on the case of magnetic fields, which is the easiest to motivate. 



t-t 



3. Magnetic field 

In the case of magnetic fields, conservation of energy- momentum, Eq. (5), has a simple 
form: e + 4de = 0. Solving this equation for d and substituting the result into Eq. (4), 
we arrive at the equation ee — + 2e^(A -|- e) = whose general solution is given 
implicitly by 

:i = i£^de(iAe^ + |efe^-e3)"^ (9) 

Once the function e is found from Eq. (9), the functions a and b are found from Eqs. (2) 
and (5) as follows: 

o=iln(q/t). (10) 

' <'^"'^ (") 

To find the behavior of a and b for large t we need to find asymptotics for the 
integral in Eq. (9). The integral diverges for small e (which corresponds to large t), 
and so we extract the divergent part first; this results in i — <j « j(3/A)2 In (ej/e) — r, 
where 

kj\e [(|Ae^)--(|Ae^-H|eJe^-e3)"'" 

Large t asymptotics give 

e~eie-^("/')*^*-*i+^\ (12) 
a- (A/3)^(i-<i + T), (13) 
b ~ (A/3) 3 (i - + r) - i In (1 + e^/X) . (14) 



^=4 
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Eqs. (13) and (14) lead to an oblate expansion with e"-^ - (1 + ej/A)2 , which agrees 
with the general result found in Eq. (7). 

We can include matter and still find exact solutions. For example, for a universe 
filled with cosmological constant, magnetic fields and dust with density p, we find 



p = Pi(e/ei)i[l + F(e/ei)]- 



A + (pi + 4ei)(6/q)^-3e 



A + pj + 



(15) 



where 



^(eAi) = |(Pi/q)[i + (A + Pi)Ai]^ 



da; a; ^ 



A/ej + (4 + Pi/ei)a;3 - 3a; 



(16) 



Equation (10) is unchanged but Eq. (11) becomes 

A+(pi + 4ei)(e/ei)t -3e 



6=iln 



A + pj + e; 

-iln(e/q) + iln[l + F(e/q)] , (17) 

and time dependence of the above functions p(e), a(e) and 6(e) can be deduced from 
the function e{t), which is now given implicitly by 



t-t: 



(18) 



4. Strings and Walls 



For strings, solutions to the Einstein equations can be obtained similarly to the 
magnetic field case. Here we give only large t behavior of these solutions: 

-2(A/3)^(/,-/.+r) 

e ~ Cj e 1 , 

a~ (A/3)^(t-ti + T), 

b - (A/3)5 (t - + r) - i In (1 + e^/X) , 



where 



T = 



de 



(iAe^)-^-(iAe^-ier^ei+e3)" 



As in the case of magnetic field, e°'~^ ~ (1 + ej/A)^, in agreement with Eq. (7). 

Solving the case of walls proceeds in a different fashion. Asymptotic forms of the 
corresponding solutions read 

-1 



-(A/3) 2 (t-t.) 

e ~ ej e i 



(l-cr)3 6aF 
1 + a ^ {1-ay 



a - (A/3)2(t -tj) - |ln(l -cr). 



6~ (A/3)5(t-fi) + ln 



(l-a)3 _^ 6aF 



1 + a (1-(t)2 
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Figure 1. The scale factors a = 6""° and P = e^~'^, where c = (|A)2 (t — t^) 
corresponds to the model with only cosmological constant. For the cases of 
magnetic field, strings, and walls, the same initial conditions axe used. Solid 
lines are exact solutions and dashed lines are asymptotics for large t. 

where 

^_ (l + q/A)^-! 
(l + ei/A)^+l 

and 




One can derive the bound F > ^{1 — a)3 {1 + a) ^, which leads to a prolate expansion 
with e"-'' < 1. 

In Fig. 1 wc compare the scale factors for the cases of magnetic field, strings, 
and walls. Since in all three cases the functions a and b asymptotically approach 
the corresponding solution for the model with only cosmological constant, c = 
(iA) 2 (t— ij), this term is conveniently factored out of the scale factors. Fig. 2 similarly 
compares the corresponding energy densities. In both figures we used the same initial 
conditions for the three cases. These conditions are fully specified by the value ej/A; 
we chose ej = A which leads to ctj = (|A)5. 

5. Conclusions 

Since the time of COBE, various speculations have been put forward to explain the 
origin of the suppressed quadrupole. One idea has been to explore topological features 
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Figure 2. The ratio e/A for the same three cases and initial conditions as in 
Fig. 1. 



with specific emphasis on the "small universe" pieture with a toroidal topology [22] . 
Another idea has been to study the possibility of a small inflation, which somehow 
gets tuned to about the scale of the present-day Hubble radius [8, 23, 24]. However, 
very little has been suggested for physical models that could explain the new features 
suggested by the WMAP data. An obvious modification of the T'^-model would be to 
make one dimension smaller than the other two and get a x asymmetric toroidal 
topology [5]. Other topologies could also work. There are many choices [26, 25] and 
the results could be profound. 

In this paper we have not departed radically from conventional cosmology, but 
rather considered a simple modification of the FRW model. We determined how 
magnetic fields, strings, and walls affect the expansion rate of the universe. Before 
solving the Einstein equations we found the asymptotic values and bounds for the ratio 
of longitudinal and transverse sizes for several interesting cases. Next, we presented 
exact analytical solutions for a universe filled with magnetic field and gave asymptotic 
results for the case of defects. Finally, we mentioned particle physics motivated 
models that could yield the sorts of magnetic fields and/or defects of importance 
for asymmetric expansion. We have found that a universe filled with aligned B-fields 
and A expands non-spherically: a co-moving sphere at evolves into an ellipsoid. 
The eccentricity is calculable and fixed by the initial conditions. Aside from being 
potentially applicable to explaining some features in the WMAP data, the models and 
the results presented in this paper differ from those in the already mentioned Ref. [12]. 
The first of the papers in Ref. [12] examined the situation where there is a uniform 
magnetic field and no cosmological constant. By adding cosmological constant, we 
have found the exact solution for the model which is physically relevant. The second 
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paper in Ref. [12] treats only small anisotropies and determines only asymptotic large 
time behavior. Our calculations apply to anisotropy of arbitrary strength and we have 
obtained the exact solutions to the Einstein equations which are valid at all times. 

The effects of non-spherical expansion and small inflation are potentially 
measurable depending on the degree of eccentricity and expansion attained and the 
quality of data. In realistic particle physics motivated models, where the B-field in 
question is associated with particular gauge fields and the cosmological constant A 
arises from some specific inflaton model, additional predictions may emerge. For 
example, a strongly aligned B-field at the beginning of infiation may affect density 
perturbations in a measurable way [27], perhaps through polarization effects, and a 
small inflation that follows may mean these perturbations are observable at the largest 
scales today. We have provided only homogeneous solutions. To apply our results to 
structure formation and CMB data, the next step will be to begin an analysis of linear 
density perturbations in such a Universe. 

The model developed here can be extended to cases of less symmetric evolution. 
Some examples are models with an aligned B-field plus cosmic strings neither aligned 
parallel nor perpendicular to this field, or models with various kinds of misaligned 
cosmic defects. Another possibility with more physical motivation is a model with 
B-fields and A where there are still a few misaligned magnetic domains within the 
horizon at the end of inflation. While this model is more complicated, it may also be 
more realistic as the initial domain size is more in line with the models of cosmological 
phase transitions [14, 15, 16, 17, 18]. Finally, while there could be other combinations 
of cosmological constant, magnetic fields, defects, and geometry that arise in branc 
world physics and produce asymmetric expansion, we believe it is important to have 
considered the simplest well motivated case first. 
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